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Abstract. - The superconducting current-phase relation (CPR) of Nb/Au/(001)YBa2Cu3Oa; 
heterojunctions prepared on epitaxial c-axis oriented YBa2Cu30i: thin films has been measured 
in a single-junction interferometer. For the first time, the second harmonic of the CPR of such 
junctions has been observed. The appearance of the second harmonic and the relative sign of 
the first and second harmonics of the CPR can be explained assuming, that the macroscopic 
pairing symmetry of our YBa2Cu30a; thin films is of the d + s type. 



Introduction. It is now well established that the dominant component of the supercon- 
ducting order parameter in the cuprates has d-wave symmetry (see Q and references therein). 
Moreover, it has become clear that in orthorhombic materials such as YBa2Cu30a; (YBCO), 
a finite component with s-wave symmetry is admixed to the dominant d-wave order parame- 
ter. The early in-plane phase sensitive experiments imply that the d-wave component remains 
coherent through the whole sample [Q, while an elegant c-axis tunneling experiment shows 
directly that the s-wave order parameter component does change sign across the twin bound- 
ary @. 

The above picture of the YBCO pairing state is challenged by the experimental observa- 
tion of a finite c-axis Josephson current between heavily twinned YBCO and a Pb counter- 
electrode jsj. Namely, the contribution of the s-wave part of the YBCO order parameter to 
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the Josephson coupling between a conventional superconductor (superconductor with a pure 
s-wave symmetry of the order parameter, for example Pb or Nb) and YBCO should average 
to zero for equal abundances of the two types of twins in YBCO. In other words, the macro- 
scopic pairing symmetry of twinned YBCO samples should be a pure c?-wave . Tanaka has 
shown that a finite second order Josephson current obtains for a junction between the s-wave 
and c-axis oriented pure d-wave superconductors However, measurements of microwave 
induced steps at multiples of hf /2e on the I-V curves of Pb/Ag/YBCO tunnel junctions 
imply dominant first order tunneling ||^. Therefore the finite c-axis Josephson current has 
to result from a nonvanishing admixture of the s-wave component to the macroscopic order 
parameter of YBCO [Q. Two alternatives of how this can take place in the junctions based 
on twinned YBCO have been discussed in the literature: 

(i) Sigrist et al. have suggested that the phase of the s-wave component in YBCO does not 
simply jump from to tt upon crossing the twin boundary, but rather changes in a smooth 
way, attaining the value of 7r/2 right at the twin boundary (t). The twinned YBCO sample is 
thus assumed to exhibit a macroscopic d + is pairing symmetry. A related picture has been 
proposed by Haslinger and Joynt, who suggest a. d + is surface state of YBCO 

(ii) A difference in the abundances of the two types of twins implies a d -I- s symmetry of 
the macroscopic pairing state Let us point out that also structural peculiarities of other 
type (such as a lamellar structure in a preferred direction) may lead to the d + s macroscopic 
pairing symmetry. 

In this paper we report the observation of a large second harmonic of the Josephson current 
in Nb/Au/(001)YBCO junctions. By comparing the relative signs of the first and second 
harmonics of the Josephson current, we show that the macroscopic d + s pairing symmetry 
is realized in our YBCO samples. Our results might be relevant also in search for so called 
"quiet" qubits which can be realized making use of junctions with a dominant second harmonic 
of the Josephson current |10|. Another promising route to fabricate such junctions is via 45° 
grain boundary Josephson junctions |Tl|p^. 

Experimental. In realization of the Josephson junctions between conventional supercon- 
ductors and c-axis oriented high-temperature ones an extensive use has been made of the 
Pb/Ag/(001)YBCO heterojunctions ||^. In search for a combination of the superconducting 
counterelcctrode and a normal-metal buffer layer with the highest possible interface quality, 
in the present work we have decided to study the Nb/Au/(001)YBCO hetcrojunction. Our 
epitaxial (OOl)-oriented YBCO films with thicknesses of 150 nm were obtained by laser de- 
position on (100) LaAlOs and (100) SrTiOa substrates. The films are usually twinned in the 
a6-plane. The superconducting transition temperature of our YBCO films was determined by 
magnetic susceptibility measurements as Tc = 88-^90 K. The YBCO films were in situ covered 
by a 8 20 nm thick Au layer, thus preventing the degradation of the YBCO surface during 
processing. Afterwards, 200 nm thick Nb counterelectrodes were deposited by DC-magnetron 
sputtering. Photolithography and low energy ion milling techniques were used to fabricate 
the Nb/Au/YBCO junctions. The interface resistance per unit area Rb — RnS (where Rn is 
the normal state resistance and S is the junction area) was Rb = 10^^ 10^^ il-cm^. Details 
of the junction fabrication were reported elsewhere ||l3|| . 

Surface quality of the YBCO films is very important when current transport in the c-axis 
direction is investigated. High-resolution atomic force microscopy reveals a smooth surface 
consisting of approximately 100 nm long islands with vertical peak-to- valley distance of w 3 nm 
(fig. 1). We can exclude that substantial afo-plane tunnel currents fiow between YBCO and Nb 
at the boundaries of these islands. In fact, theory predicts formation of midgap states at the 
surface of semi-infinite Cu02 planes ||lj, Therefore zero bias conductance peaks should 
be expected in the I-V characteristics at temperatures larger than the critical temperature 
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of Nb, if the contribution of ab-plane tunneling was nonnegligible. However, no such peaks 
have been observed for all fabricated Nb/Au/YBCO junctions. Moreover, from the size of the 
islands and from the vertical peak-to-valley distance we estimate that the area across which 
a5-plane tunneling might take place is less than 6% of the total junction area. Since the 
interface resistances per square are of the same order of magnitude ^ for both, the c-axis and 
a6-plane junctions, we conclude that a6-plane tunneling from YBCO, if present, is negligibly 
small. 

We have measured more than 20 single junctions with areas in the range from 10 x 10 /im^ 
to 100 X 100 jjLw? . At small voltages the typical I-V curves can be described by the resistively 
shunted junction model with a small capacitance p6[ |. Typical critical current densities were 
jc = 1-^ 12 A/cm^ and jcRs — lO-f-OO /iV. The differential resistance vs. voltage dependence 
Rd{V) exhibits a gap-like structure aXV~ 1.2 mV at T = 4.2 K (see fig. 2). This structure 
has a BCS-like temperature dependence and disappears at Tcr ~ 9.1 K, therefore we ascribe 
it to the superconducting energy gap of Nb. 

The CPR measurements were performed using a single-junction interferometer configura- 
tion in which a junction of interest is inserted into a superconducting loop with an inductance 
L « 80 pH. We measure the impedance of a parallel resonance circuit inductively coupled to 
the interferometer as a function of the external magnetic flux <l>e threading the interferome- 
ter. The dimensionless CPR fi^p) = Is{<f)/Ic (where Is{^) is the Josephson current) can be 
extracted from the following equations: 

^ = V'e-ZJ/M, (1) 

= YTm^r 

where (p^ = 27r$e/<I'o, is the phase difference across the Josephson junction, the prime 
denotes a derivative with respect to if, a is the phase shift between the driving current and 
the tank voltage at the resonant frequency, (3 = 2nLIc/^o is a normalized critical current, Q 
is the quality factor of the parallel resonance circuit, k is the coupling coefficient between the 
RF SQUID and the tank coil, and $o is the flux quantum. This method, being differential 
with respect to (p, provides a high sensitivity of the CPR measurement. Moreover, Is{ip) is 
measurable even if the thermal energy exceeds the Josephson coupling energy. In fact, critical 
currents down to 50 nA were recently detected at T = 4.2 K pTf . 

For the junction SQIO, /? varied between 0.4 and 0.27 for temperatures in the range 
T = 1.7 6.0 K. Since (3 < 1, we can extract the CPR from the a(</5e) dependence for the 
complete phase range (see fig. 3). Fourier analysis of the experimentally obtained CPR shows 
substantial first and second harmonics and negligibly small higher-order harmonics. Therefore 
we can write 

Is ip) = h sin p + h sin 2p (3) 

for all temperatures below the transition temperature of Nb. The sign of I2 is always opposite 
to that of Ii. In what follows we use the convention that Ii > and therefore /2 < 0. The 
ratio l/2/^il grows with decreasing temperature reaching I/2/AI ~ 0.16 at T w 1.7 K, when 
Ii = 1.57 fiA and I2 = —0.25 ^A (see inset to a fig. 3(b)). 

We point out that the main result of this paper, namely the large negative value /2//1 ~ 
—0.16 observed at low temperatures, is not the result of an indirect data analysis, but it 
follows directly from the measured a{p>^) dependences. In fact, one finds readily that 

rf^tang _ k'^QPf" 
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where all derivatives are taken a,t tp — or ip^^ — 0. Thus, the existence of the local minima 
of at multiples of tt dictates that /"' > 0, or, making use of Eq. (|), h/h < -1/8. 

Note that neither the conventional tunneling theory nor the I and II theories by Kulik and 
Omelyanchuk predict such local minima on the derivatives of CPR [ p^ . 

Discussion. Let us estimate the transparency of the barrier between YBCO and Nb from 
the normal-state resistance per unit area Rb- According to the band-structure calculations 
(for a review, see the hole Fermi surface of YBCO is a slightly warped barrel with an 

approximately circular in-plane cross-section (to be called Fermi line) with radius kp. In what 
follows, we represent the electron wavevector k in cylindrical coordinates, k = (fc, 0, k^). We 
estimate the uncertainty of the in-plane momentum as 5k w 27r//, where I is the characteristic 
size of the islands on the YBCO surface (fig.l). We evaluate Rb making use of the Landauer 
formula and note that only tunneling from a shell around the Fermi line with width Sk is 
kinematically allowed. The barrier transparency D{9) depends on the details of the c-axis 
charge dynamics in YBCO, with maxima in those directions 9, in which the YBCO c-axis 
Fermi velocity w{6) is maximal. Since for 6 = 7r/4 and symmetry equivalent directions w{9) is 
minimal [|l9|, we expect that there will be 8 maxima of D{9) on the YBCO Fermi line where 
D(9) ~ D, which are situated at 9 = 9o and symmetry equivalent directions. The modulation 
of the function D{9) depends on the thickness of the barrier between YBCO and Nb We 
consider two limiting distributions of the barrier transparency D{9) along the YBCO Fermi 
line: (a) a featureless D{9) w D and (b) a strongly peaked D{9), roughly corresponding to 
thin and thick barriers, respectively pO| ]. In the thick barrier limit the angular size of the 
maxima of D{9) can be estimated as 69 ~ 6k/ kp. With these assumptions we find 



where A measures the number of conduction channels and (. . .) denotes an average over the 
junction area. In the thin and thick barrier limits, we find A « kpSk/n and A « 26k'^ /tt, 
respectively. Taking I k. 100 nm and kp ~ 0.6 [^, the measured Rb = 6 x 10~^ Vtcw? 
of the junction SQIO can be fitted with {D)thin ~ 1-7 x 10"^ and (i:')thick ~ 8.3 x lO"'^. 

Since we have observed no midgap surface states in the Rd{y) curves, we can neglect the 
surface roughness, and the Josephson current can be calculated from ||2^ 

7 foo^ = ^ Va- r V UAj7. Ak sin Lp 

^ It ^ 4;^ 2f]fl,l]k + D[lo'^ + VLrQ.^ + Afl. Ak cos ifY ^ ' 

where the sum over k,9 is taken over the same regions with areas A as in Eq. (||), A/j and 
Ak are the Nb and YBCO gaps, respectively, and fli — y/uj'^ + Af with i = R,k. Keeping 
only terms up to second order in the (small) junction transparency D, the Josephson current 
densities ji = li/S read 

3i{T)Rb ~ -TT , (7) 

32{T)Rb ~ ~^-rFT\ : — t^^^ ni. rr. , (8) 



8 (D) e V 2fcBT 

where AjJ = 7rA(;[2 ln(3.56A(;/rcfl)]^^ and A^ = Ad|cos26'o| in the thin and thick barrier 
limits, respectively. In Eqs. (0,||) we used the YBCO gap A{9) — As + A^ cos 20, where 
Ad and A^ are the d-wave and s-wave gaps. We have assumed that A^ is larger than both. 
An and A^. The factor As/A^ can be estimated from the measured JiRb products for 
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Josephson junctions between untwinned YBCO single crystals and Pb counterelectrodes. For 
such junctions ji(0)i?s ~ 0.5 1.6 meV Using the Pb gap A^j = 1.4 meV in Eq. (|^), we 
obtain A^/A^ ~ 0.36 1.1. 

Note that within the above microscopic d + s scenario, the signs of Ii and I2 are different. 
This feature remains valid also in the macroscopic d + s scenario, whereas within the macro- 
scopic d + is picture, the same signs of /i and I2 are expected. Thus we conclude that the 
finite first harmonic has to be due to the macroscopic d -\- s symmetry of our YBCO sample. 
In fact, detailed structural studies show that for sufficiently thin YBCO films, the abun- 
dances of the two types of twins can be different even for films grown on the cubic substrate 
SrTiOs |2^. If we denote the twin fractions as (1 -I- 6)/2 and (1 — (5)/2, then the measured 
first harmonic of the CPR, (ji), is proportional to the deviation from equal population of 
twins, (ji) = (5ji Using A^j = 1.2 meV determined from the Rd{V) data and our estimate 
Ag/AjJ K, 0.36-^ 1.1, we find that the measured first harmonic (ji) for the junction SQIO can 
be fitted with 5 « 0.07 0.21, which is in qualitative agreement with where 5 ~ 0.14 for 
1000 A thick YBCO films has been observed. 

Fitting the measured j2-R_B of the junction SQIO by Eq. (||) we obtain (D^)/(_D) « 3.2 x 
10"^, which is much larger than both (I?)thin and (£')thick- This diference can be explained 
provided the junction transparency _D is a fluctuating function of the position r. In fact, 
adopting the WKB description of tunneling ||2^, we write D{s{y)) — exp(— sq — s(r)), where 
So is the WKB tunneling exponent and s(r) its local deviation from the mean. Assuming a 
Gaussian distribution of s with a mean deviation r/, P{s) cx exp(— s^/77^), we estimate the 
spatial averages as (-D") = /^°^ dsP{s)D'"-{s). In the thin barrier limit, the values (£'^)thin = 
8.6 X 10^^ and (£')thin = 2.3 x 10^^ required to fit the experiments correspond to an average 
WKB exponent Sg'''" « 15.5 with jythin ~ 4-3. In the thick barrier limit we obtain Sq^"^'*^ k, 9.1 
and 77thick ~ 2.8. 

Let us consider also the second harmonic generation by a mechanism proposed by Millis for 
planar junctions p^ . We can view the junction as a checkerboard of and tt junctions with a 
lattice constant a (characteristic size of an YBCO twin) and local critical current density ji. 
As shown by Millis, spontaneous currents are generated in the ground state of such a junction 
and the junction energy is minimized for the phase difference ±7r/2. An explicit calculation 
in the limit a, A^j ^ A ^ Ac (where A, Ac are the a5-plane and c-axis penetration depths of 
YBCO and Xr is the Nb penetration depth) yields 

1 fi^jfaXXc 

J2, Millis 



4\/2 $0 

Comparing j2,Miiiis with Eq. (||), we find 

j2,Miiiis 8V2F f As\^ a Ar Eo 



32 7r2 \AlJ Ific/Xcfic/y 

where Eq = e^qo/Aireo with 50 — kp and 26k for a thin and thick barrier, respectively. F 
depends on the characteristic length scale L of the fiuctuations of the junction transparency: 
ioi L <^ a, F Ki {D)'^/{D'^), whereas for L ^ a, F « 1. In order to estimate the upper 
bound of the ratio j2.Miiiis/j2, we consider the thin barrier limit, take F Ri 1, As/A^ « 1.1, 
a « 10 nm, A ~ 240 nm, and Ac ~ 3 /xm (the last two values are valid for underdoped 
YBCO [^) and we find that j2,Miiiis/i2 < 0.03. Thus we can safely neglect the contribution 
of the Millis mechanism to the second harmonic of the CPR. 

Conclusion. We have observed the second harmonic I2 of the current-phase relation in c- 
axis Nb/Au/YBCO heterojunctions. We have shown that the relative phases of the first and 
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second harmonics of the CPR together with their amphtudes Ii and I2 and the normal-state 
resistance of the Josephson junctions observed in the experiment can be explained making 
use of a single set of parameters within the standard microscopic d + s picture of the YBCO 
pairing state, assuming the difference in the abudances of the two types of twins and that the 
barrier fluctuates along the junction. 
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Figure captions 

Fig. 1. Left panel: A high resolution AFM image of the surface of a 150 nm thick YBCO 
film. Right panel: Height profile of the YBCO surface along the line indicated in the left 
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panel. The peaks and valleys are shown by markers. The peak-to-valley distance is ~3 nm 
and about 100 nm in the vertical and horisontal directions respectively. 

Fig. 2(a) Typical I — V curve and (b) differential resistance vs. voltage dependence Rd{V) 
of a Nb/Au/YBCO junction with a 20 nm thick Au film, measured at T = 4.2 K. The low 
voltage part of the I — V curve ( V < 0.3 mV) is shown on the inset. 

Fig. 3. Phase shift a as a function of ip^ for the jimction SQIO at T =1.7, 2.5, 3.5, 4.2, 
and 6.0 K (from bottom to top). 

Fig. 4. The current-phase relation I{(p) of the junction SQIO at T=1.7, 2.5, 3.5, 4.2, and 
6.0 K (from top to bottom). Inset: Temperature dependence of Ii (squares) and I/2I (circles). 
Solid lines are fits to Eqs. (|,|) using Afl(T) = Aij(O) tanh[Aij(r)Tcfl/Afl;(0)T]. 




Fig. 4. 



